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Bohlin- Arnold- Vassiliev's duality transformation establishes a correspondence between motions 
in different central potentials, ft offers a very direct way to construct the dynamical conserved 
quantities associated to the isotropic harmonic oscillator (Fradkin-Jauch-Hill tensor) and to the 
Kepler's problem (Laplace-Runge-Lenz vector). 
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(N ■ I. INTRODUCTION 
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The study of the relation between the classical motions of the harmonic oscillator and the Kepler problem has a long 
history initiated in Hooke and Newton's worksA. In 1911, Bohlin showed that this relation can be formulated in terms 
of conformal mapping^. This conformal transform is at the core of the Levi-Civita's regularization scheme 3 - of which 
tri-dimensional generalization has been achieved by Kustaanheimo and Sticfcl 4 . In fact, as noted by Needham££ two 
years before Bohlin's paper, Kasner— had established a more general duality law relating pairs of power law potentials 
(but the result has only been published in 1913). This relation has been rediscovered and generalized by Arnold and 
'-"^ ■ Vassiliev^ in 1989 and quite simultaneously by Hojman and alJ£. In fact, the generalization of Kasner's result had 
II | already been obtained by CdHas 1 *^ and implicitly enters into the frame of the coupling constant metamorphosis of 
Hietarinta and a h 12 ' 13 i 14 . Even if we limit ourselves to the classical aspects (The study of this correspondence in the 
quantum mechanical frame has an interesting parallel history that we won't deal with here), numerous articles have 
been published on this subject during the last fifteen year o 1 ! 15 ! 16 ! 17 ! 18 ! 19 ! 20 . In all what follows, we stay in a purely 

1 i ' Newtonian frame. On the basis of simple but general arguments concerning the complex representation of 2D-motions, 

we show that the Bohlin- Arnold- Vassiliev's dual correspondence appears very naturally as the relevant transformation 
O^l ' of motions among all combinations of analytic change of coordinates and Euler-Sundman reparameterization. We 
consider the effect of this duality on conserved quantities. In the specific cases of Hooke's and Kepler's problems, 
we recover in a very simple and direct way the correspondence established by Nersessian and ah 28 ' 29 i 30 between the 
associated additional dynamical conserved quantities. This permits to point out the essentially obvious character of 
. Fradkin-Jauch-Hill tensor's and Laplace-Runge-Lenz vector's conservation. 



II. MOTIONS TRANSFORMATIONS 
A. Complex formulation 
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Let's consider a planar motion V (t) = (^n) _» _ for a particle of mass m submitted to a potential iK^r), 

eventually singular at the origin. We will adopt a complex formulation and represent the position by its corresponding 
affix z(t) = x(t) + iy(t), the potential being viewed as a real valued function, t/(z,z). The gradient of such a real 
valued function of z and z, is given by: 

w ( y) = 2^M a) 

where ■§= = \ (j^ + i^j = ~f z . In particular, for a central potential: U(z,z) = U (\z\) = U (r) and = ^p-f. 

The equation of motion takes the form: 

. + A^I^O (2) 
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where the dot represents the time t derivative. 

Multiplying both sides by z and taking the real part of the resulting quantity, we obtain a total derivative: 
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Integrating this identity we recover the conservation of energy: 
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If yl and i? are two vectors in the (O, u^, uj) plane, the complex affix of A x i3 is given by the real quantity Im (AB^ , 

A and B being the complex affixes of A and B respectively (Note that this real quantity doesn't correspond to a 
vector in the (0,u x ,u y ) plane but which is orthogonal to this last). 

The angular momentum L (t) = m r (t) x r (t) = L (t) ul admits then as complex correspondent L (t) 

in [in (z(t)z{t)) =f i (zz- zz^j . 

Moreover: 



L{t) = Y l ( zz ' zz ) 

which implies that for a central potential L (t) = 0. 



f_dU(z,z) dU(z,z) 
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B. Conformal change of coordinates 



Consider first an arbitrary conformal change of coordinates: 



Then 



= f(w), z = f(w) = f(w) 
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where f^(z) = 

After substitution, Eq.([21) takes the following form: 



(6) 
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w + (w) 2 log (/ (1) H) 
where U(w,w) = U (f(w),f(w)). 



(!) 2 1 dU(w,w) 



in 



\fW(w)\ 2 dw 



(8) 



C. Euler-Sundman's reparameterization 

The previous equation {5| incorporates now a term of first order in time derivative, which has an artificial character. 
In order to suppress it, we can envisage to change the parameterization for the motion w(t). A global change of 
parameterization will be inefficient. Rather, we choose a change of parameterization having a local structure, that is 
an Euler-Sundman's reparameterization of the type: 



dt=\g(w)\ 2 ds (9) 

where g (w) is an, a priori, arbitrary analytic function. This form ensures in particular that the correspondence 
between the initial and the new (fictitious) time is one-to-one and increasing. 
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Starting from Eq.®, it is straightforward to show that the equation for the motion w(s) is now given by: 

W " + K) 2 (log ^) W - (io S9 H)^ W\ 2 + - ^L 9 ^ 1 = o do) 

V g(w) J m |/(i)( w )| ow 

A priori, the structure of this last equation is still more complicated than the previous one, since it incorporates 
two kinds of terms depending on first derivative in time: the first one depending on (w 1 ) 2 as in Eq.®) and the second 
one, depending on \w'\ . But the gain becomes readily clear. Indeed, the freedom in the choice of g (w) allows to 
expect to suppress the first one. As for the problem induced by the second, as we just go to see, it is solved by the 
use of energy conservation identity. 



D. Energy conservation 

Expressed in terms of w and s, the energy conservation equation gives: 



i„f = 2( E -5 ( „,w))-M^ an 

/' ■ >(w)\ 
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We are then in position to replace in Eq. (|10p the contribution containing \w'\ 2 by another one, depending only of 
the position. We obtain: 

E. Reparameterization's choice 

In order to eliminate the contribution containing (w') 2 , we have to fix g (w) such that g (w) = C/^ 1 ' (w) , C G C. 
The simplest choice, which we'll adopt in the sequel, is to take C = 1 and: 

g(w)=f^(w) (13) 

Then Eq.(??) becomes simply : 

w" + ^f [1 \w)^(j^)(u(w,w)-E)) =0 (14) 

Defining : 

V f (w,w)= f (1) (w) 2 (u(w,w) -E) +V (15) 
where Vo is an arbitrary real constant, we finally obtain : 

w „ + 2, dV f (w,W) =Q 
m ow 

In other words, w (s) is the motion of a particle of mass m in the potential Vf(w,w) (if we choose as supplementary 
additive term an arbitrary holomorphic function Vo(w), we are led to the same equation but the corresponding 
potential Vf(w,w) is no more real valued). 

As for the energy of this dual system, it's given by : 

E= - \w'\ 2 + V f (w,w)=V (17) 
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Relation (fl~5|) becomes finally: 

(1) 2 = E-Vj(w,w) = E-V f (w,w) 

E-U(w,w) E-U(f(w),f(w)) 

We then obtain a functional equation to which is submitted the transformation / linking the motions in the 
potentials U and V respectively. 

In conclusion, among all the transformations of a 2D-motion in a potential U combining a conformal (i.e. analytical) 
change of coordinates z — f(w) and an Euler-Sundman reparameterization dt = \g (w)\ 2 ds, those for which the 
transformed motion is an autonomous 2D-motion submitted to a potential V, necessarily satisfy Eq. fHi]) and Eq. (fT5)l . 



III. BOHLIN-ARNOLD-VASSILIEV'S DUALITY 



A. Arnold- Vassiliev's potentials 



We consider the case where the initial potential Ui^r) takes the form Ui^r) — U(z,z) — Au(z)u(z) — A \u(z)Y , 
where u(z) is an analytical function on C* and Ael. We'll call such a potential, an Arnold- Vassiliev's potential 9 . 
It is easy to show that the only central Arnold- Vassiliev's potentials are the power law potentials, corresponding to 
u(z) ~ z% , that is: 



U(z,z) = A\z\ v , i/,4el (19) 

If u(w) = (u o /) (w) , we have, by the transformation associated to /, an associated dual potential of the form (see 
Eq.©): 



E-V{w,W)= (E-A\u(w)\ 2 ^) 



(20) 



Let's ask for which type of conformal transformation /, the corresponding image potential V is also an Arnold- 
Vassiliev's potential V(w, w) — B \v(w)\ . 

Inserting in Eq. (|2"0)) and deriving successively with respect to w and w, we obtain: 



AE\v{w)\ 2 (E-B\v{w)f 



BE \u(w)Y 



E-A\u{w)Y 



(21) 



If we choose B and E such that = 1, this gives the following simple constraint : 



\u{w)\ 



/ (1) H 



(22) 



which is clearly satisfied when: 



v(w) = fU(w) = -r—=— (23) 
u(f(w)) u(z) 



that is when w — f 1 (z) = J u(z)dz. 
Then Eq.® gives: 



E-B\fW(w)\ a 
E\fW(w)\ 2 -A 



(24) 
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which, jointly with 4tI = 1, implies: 

BE 



E = -A,B = -E (25) 

The transformations satisfying Eq. (|23|) and Eq. (|25|) let stable the set of Arnold- Vassiliev's potentials. Since they 
are involutions, they can be seen as duality transformations: 



(E,A\u(z)f) W =V d * ("4 j^jj) ( E ,A\u(z)f) (26) 

The Arnold- Vassiliev's potentials appear then as the more natural class of potentials for which the above combination 
of an analytical change of coordinates and an Euler-Sundman reparameterization is a dual correspondence. 



B. Power- law potentials 



Let's consider the particular case of power-law potentials: u(z) = z^, U(z,z) = k \z\ , at fixed energy E. Eq.f2|) 
becomes: 



— \z\"-*z = (27) 
m 



The preceding results (see Eq. (|26|0 give: 



w = f-\z) = T ^z^\ z = f{w)=(l + ^)^ (28) 

By the duality transformation considered above, the corresponding image potential (see Eq.(??), Eg. ([25]) and note 
that A = k) is then given by: 



v(w) = (l + ~) '" +2 w-^, V(w,w) = B\v(w)\ 2 = k\w\" (29) 

where: 



" = -if|^=- E ( 1 + f)" < 3 »> 

The energy for the dual motion is (see Eq. (|25]) ) E = — fc and the equation of the dual image motion takes the form: 



w 



— \w\ w 

m 



(31) 



We recover for the dual image motion, a motion of energy (—A), submitted to a power law potential of characteristic 
exponent /i, such that: 



( 1+ D( 1+ fH ^ 

and for which the coupling constant is proportional to the opposite of the initial energy E (and real if v > —2). 
The relation between the two motions will be called "Bohlin- Arnold- Vassiliev's duality". 

In the special case v = 2, we obtain fj, = — 1. In other words, the dual motion of the plane harmonic oscillator 
(Hooke potential ^kr 2 , energy E) is nothing but the Kepler motion (Newton potential — energy — §)• We recover 
the usual Levi-Civita regularizing transformation with the associated change of variables: 
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z = \/2w, w = —z 2 (33) 

and reparameterization: 

ds = 2pdt = r 2 dt (34) 

where r = \z\ and p = \w\. 

IV. DUALITY AND CONSERVED QUANTITIES 

For a general 2D-central potential we have two real conserved quantities: the energy E, associated to translational 
invariance in time, and the angular momentum L, associated to rotational invariance. In the case of power-law poten- 
tials, the effect of the duality transformation on the energy is given by Eq. ([3H| . Concerning the angular momentum, 
it is easy to show that (see Eq. (p~3|) ) : 

171 / . j_ \ Tn ( w' w' 



n (_. _\ 

— zz — zz ] 



2iK ' 21 V(io g /H) (1) (iog/H) (1 V 

which gives for power law potentials (see Eq . (f58|) ) : 
where L is the angular momentum of the w{s) motion. 

Among all the central potentials, only the Kepler system and the isotropic harmonic oscillator possess true additional 
conserved quantities which are the Laplace-Runge-Lenz vector— and the Fradkin-Jauch-Hill tenso r 22 ! 23 respectively. 
With all the other central potentials are associated only piecewise conserved vector o 24 ' 25 . 

This peculiarity of Kepler and harmonic potentials is at the origin of the existence of closed orbits for all the bound 
state a 26 i 27 . 

In complex representation the existence of a specific conserved quantity for the harmonic oscillator becomes a 
matter of course. Indeed, if we consider the equation for a plane motion in a central potential U(r): 

1 U {1 Hr) 

z + -- ^z = (35) 

m r 

We see immediately that among all these potentials, the only one for which the corresponding equation is not z 
dependent is the linear one, that is the one associated to the Hooke potential (i.e. isotropic harmonic oscillator): 

k 

z+—z = (36) 
m 

This characteristic feature allows the existence of an immediate integrating factor z, which after integration leads 
to the following complex conserved quantity: 

/ ■ \ 2 2 

T = mA^— + k— = cste (37) 

We then have: 



T = (Tn - T 22 ) + iT 12 



(38) 
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where Ta is the well known Fradkin-Jauch-Hill's tenso r 22 i 23 : 

Ti | TV, \ _ / f f.rr + fz 2 mf + fcf 



l V7 

v 2 



As for the energy, it is nothing else but the trace of this tensor: 



(39) 



E = T n +T 22 = ^m\z\ 2 + ^k\z\ 2 (40) 

In complex representation, the existence of Fradkin-Jauch-Hill's tensor appears therefore as a trivial consequence 
of the linearity of the equation of motion. 

Let's interest ourselves to the interpretation of T for the associated dual Kepler motion. In terms of w variable and 
s parameter (see Eq. ([33l) and Eq. (f34|) ). T takes the form : 

T = — — (w) 2 + kw = rriw(w') 2 + kw (41) 

Since the energy of the dual motion is given by E = — -| = |u/| 2 — we can write: 

T = mw' (ww - wvf) + E— = 2imw'L - 2k— (42) 

P P 

where k = — 

If we note p = £/u x + rfu y the vector of affix w, iLw' is the affix of L x p ', We then obtain: 

T = 2kA = -EA (43) 

where: 

A = —iwL (44) 

k P 

We immediately recognize the affix of the well-known Laplace-Runge-Lenz vector— for the Kepler motion w: 

-Z = ™Z x y>_l (45 ) 

k ' P 

By Eq. (|43| we then have a very direct interpretation of the Laplace-Runge-Lenz vector as the ratio of two conserved 
quantities associated to the dual motion of isotropic harmonic oscillator: 

A = -\ (46) 
This ensures immediately the conservative character of A. 

The complex Newtonian frame adopted here permits to bring out the obvious character of the conservation Laplace- 
Runge-Lenz vector as a consequence, via the dual transform, of the triviality of the Fradkin-Jauch-Hill tensor's 
conservation. It is interesting to note that the elementary approach developed here can be extended to a much more 
general set of planar motions conserving only the angular momentum directional. In this case, the duality connects 
different classes of motions. For those presenting closed orbits which satisfy generalized Gorringe-Leach equations, 
there exist additional conserved quantities which are linked by a relation similar to Eq. 
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